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Abstract 

A macroscopic model of the tumor Gompertzian growth is proposed. This approach is based 
on the energetic balance among the different cell activities, described by methods of statistical 
mechanics and related to the growth inhibitor factors. The model is successfully applied to the 
multicellular tumor spheroid data. 



I. INTRODUCTION 



A microscopic model of tumor growth in vivo is still an open problem. However, in spite 
of the large set of potential parameters due to the variety of the in situ conditions, tumors 
have a peculiar growth pattern that is generally described by a Gompertzian curve often 
considered as a pure phenomenological fit of the data. More precisely there is an initial ex- 
ponential growth (until 1-3 mm in diameter) followed by the vascular Gompertzian phase Q]- 
Then it seems reasonable to think that cancer growth follows a general pattern that one can 
hope to describe by macroscopic variables and following this line of research, for example, the 
universal model proposed in 31 has been recently applied to cancer j^. In this talk we present 
a macroscopic model of tumor growth, proposed in ^, that: i) gives an energetic basis to 
the Gompertzian law; ii) clearly distinguishes between the general evolution patterns, which 
include the internal feedback effects, and the external constraints; iii) can give indications 
on the different tumor phases during its evolution. The proposed macroscopic approach is 
not in competition with microscopic models j^, but it is a complementary instrument for 
the description of the tumor growth. 



II. CELLULAR ENERGETIC BALANCE AND GOMPERTZIAN GROWTH 

The Gompertzian curve is solution of the equation 

where N{t) is the cell number at time t, 7 is a constant and N^ax is the theoretical saturation 
value for t ^ 00. 

It is quite natural to identify the right hand side of Eq. (Q) as the number of proliferating 
cells at time t and then to consider fp{N) = 7 In ^'j^"" ^ as the fraction of proliferating 
cells and 1 — fp{N) = fnp the fraction of non proliferating cells. Since fp{N) depends on 
N{t), there is a feedback mechanism usually described by introducing some growth inhibitor 
factors which increase with the number of non-proliferating cells and are responsible for the 
saturation of the tumor size. The concentration of inhibitor factors should be proportional 
to the number of non-proliferating cells which is maximum at = N^ax Q| • If one considers 
that, during the growth, each cell shares out its available energy at time t, in the average, 
among its metabolic activities, the mechanical work ( associated with the change of the 
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tumor size and shape) and the increase of the number of cells, it is conceivable to translate 
the previous cellular feedback effect in terms of energy content. Indeed, as shown in P], the 
specific energy for the growth should be proportional to fnp and the average metabolic plus 
mechanical energy per cell, Mg, is proportional to fp. As we shall see this reproduces the 
observed cellular feedback. 

The model P], based on an analogy with statistical mechanics, assumes that in a larger 
system B, the body, there is a subsystem A, the tumor, made of N{t) cells at time t, with 
total energy U, which has specific distributive mechanisms for providing, in the average, 
the amount of energy U/N to each cell. Then we indicate with Em the energy needed 
for the metabolic activities of A, with Q the energy associated with the mechanical work 
required for any change of size and shape of A, and with /i the specific energy (i.e. per 
cell) correlated to the change in the number of cells A^, and, by assuming that these three 
processes summarize the whole cellular activity, we have U = Em + ^ + ^iN . Let us assume 
that the system A slowly evolves through states of equilibrium with the system B, defined 
by macroscopic variables, analogous for instance to the inverse temperature /5, that have the 
same value for the two systems, although it should be clear that in our case we do not have 
real thermodynamical equilibrium because the system B supplies the global energy for the 
slow evolution of the subsystem A. Within this scheme, there are many microscopic states 
of the system A, compatible with the macroscopic state, defined by /3, ^ and V , which are 
built by a large number of states of each single cell. These microscopic states of each cell 
have minimum total energy e and in an extremely simplified picture, we assume an energy 
spectrum of the form = e + 15, where / is an integer and 5 is the minimum energy gap 
between two states. With this spectrum the grand partition function Z is given by the 
following product Z{(3,V,iJ,) = II^Qexp ^^e"^*^*^'"^^^ and the corresponding grand potential, 
which is natural to associate to the energy Q related to the mechanical work in our problem, 
is given by fl{l3, V, /i) = — (i?//9)exp[— /5(e — fi)] where R= 1/(1 — e~^^). The average value 
of A^, defined for constant V and jS, turns out to be A^ = —l3Q. According to the basic rules 
of statistical mechanics, the product of the "entropy times the temperature", which in our 
system corresponds to Em introduced above, is 



where C is given by C = RP6 exp(— /5(5). From the previous equations it is straightforward 
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to express /i in terms of : /i = e + ln(A^/_R). 

To find the evolution of the system A which takes into account the internal feedback 
mechanism we recall some results obtained in Isf : 



1. the energetic balance requires that the growth with cellular feedback starts at a 
minimum number of cells Nm and saturates at a maximum value, N^ax related by 

Nmax = iVmexp(l + pe). 

2. For Nm » 1, Me = Em/N + Q/N = (1/(3) \n{Nmax/N) is a decreasing function of 

> Nm and there is a simultaneous reduction of the total metabolic energy per cell 
and an increase of specific energy required for the growth: there is an energetic balance 
between Mg and fi = (1//5)(1 + /?e + In^N / Nmax)) ■ Moreover Mg is proportional to 
fpit). 

Then it is possible to derive the Gompertz equation for the growth, AA^ in an interval 
At. AA^ is proportional to the number of proliferating cells and then one can write AA^ = 
Ci At fp{t)N = C2 At Me A^ where ci and C2 are constants. This gives in the continuum 
limit Eq. (P) with 7 = C2/P. 

III. APPLICATION TO MULTICELLULAR TUMOR SPHEROIDS 

The first step to analyze the phenomenological implications of the model and to describe 
the dependence of the growth on the external conditions is to consider the multicellular 
tumor spheroids (MTS). A minimal MTS description consists of a spherical growth where 
a) the thickness, k, of the layer where the nutrient and oxygen is delivered (the crust) is 
independent on the spheroid radius R; b) the cell density is constant; c) the cells in the crust, 
receive a constant supply of nutrient for cell; d) at time t the cells are non proliferating if 
they are at distances d<R — k if R>k from the center of the spheroid. For R < k all cells 
are proliferating. 

To separate the effects of the external constraints due to energy supply from those related 
to biomechanical conditions, it is better to consider first the MTS growth without external 
and internal stress and to introduce later these effects. 
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A. Energetic MTS growth 



In this case the external conditions are experimentally modified by changing the oxygen 
and nutrient concentration in the environment. At fixed value of these concentrations, the 
maximum allowed number of cells in the MTS is N^ax- For R < k dX\ cells receive the 
nutrient and oxygen supply while for R > k there is a fraction of non proliferating cells and 
the feedback effect starts. The growth of the MTS is due to the proliferating cells in the 
crust and one obtains that the MTS radius R, for R >> k, follows a Gompertzian law as 
the one in Eq. (^J), with replaced by R and N^ax by Rmax, where Rmax is the maximum 
radius of the spheroid corresponding to the maximum number of tumor cells N^ax- The 
experimental results show that after 3-4 days of initial exponential growth the spheroids 
essentially follow the Gompertzian pattern 7|. 

According to our model, at time t > t*, such that R{t*) = k, a fraction of the total cells 
becomes non proliferating, the feedback effect starts and the growth rate decreases according 
to the Gompertz law. The number of cells at time t* is fixed by the condition N{t*) = Nm p|- 
On the other hand, the variation of the concentration of nutrient and/or of oxygen modifies 
the total energy supply, that is the value of N^ax and, since = Nmax^w{~^P ~ 2), 
there is a clear correlation among the external energetic "boundary conditions", the value 
Nmax and the thickness of the viable cell rim which corresponds to the radius of the onset 
of necrosis. It can be shown j^] that {Gc is the glucose concentration): 

k{G,) = a {n'JI - N'^J^') + ko (3) 

where a and k^ are constants depending on the supplied oxygen. From Eq. one obtains 
the correlation among N^ax, Gc and k. In Fig.l and Fig. 2 the previous behaviors are 
compared with data without optimization of the parameters. 



B. Biomechanical effects 



The experimental data indicate that when MTS are under a solid stress, obtained for 
example by a gel, the cellular density p is not constant and depends on the external gel 
concentration Cg. In particular the results in show that: 1) an increase of the gel 
concentration inhibits the growth of MTS; 2) the cellular density at saturation increases with 
the gel concentration. In the model the mechanical energy is included in the energetic balance 
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Cg (percent) 2Rmax{^) [fJ-m] exper. 2Rmax{^) [f^m] fit 



0.3 



450 



452 



0.5 



414 



429 



0.7 



370 



404 



0.8 



363 



394 



TABLE I: Comparison with the experimental data as discussed in the text. The experimental error 
is about ±10%. 

of the system by the term Q = ~N/f3 = —PV where the pressure is P{t) = p{t)/f3. The 
introduction of this term decreases the value of Nmax with respect to the case in Sect. IIIlXl 
and this reduction should also imply a decrease of the maximum size of the spheroids, i.e. 
Rmax{P) by increasing the pressure. The comparison with the data is reported in Table I 
for Cg in the range 0.3 — 0.8 % (see i^.] for details). 
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FIG. 2: Spheroid saturation cells number vs. diameter (/xm) at which necrosis first develops. 



Circles refer to culture in 0.28 mM of oxygen Triangles refer to culture in 0.07 mM of oxygen. 
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